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Abstract 

We demonstrate how the virtual photoproduction of vector mesons on nuclei scans 
the wave function of vector mesons from the large non-perturbative transverse size p ~ 
Rv down to the small perturbative size p ~ l/^/Q^. Thee mechanism of scanning is 
based on color transparency and QCD predicted spatial wave function of quark- ant iquark 
fluctuations of virtual photons. A rich, energy- and Q^-dependent, pattern of the nuclear 
shadowing and antishadowing is predicted, which can be tested at the European Electron 
Facility and SLAC. 
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Introduction. In this paper we discuss the novel feature of color transparency (CT) 
tests of QCD: the scanning of the nonperturbative hadronic wave functions. The virtual 
photoproduction of vector mesons (T, T', J/^, p, p', ...) is particularly suited for such 
CT tests. At high energy z/, the photoproduction can be viewed [1,2] as a production of 
the virtual qq pair with the coherence length 
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Ic^ , o • (1) 



The size pq of the produced qq pair (the ejectile state) is controlled by the virtuality of 
photons [3]. The qq state is projected onto (recombines into) the final-state vector meson 
V with the formation (recombination) length 

where Am is the typical level splitting in the quarkonium. We demonstrate, how by 
changing and the initial size pg one can scan the wave function of vector mesons 
starting from the non-perturbative size p ~ Ry down to the perturbative region of p ~ 
1/Q. 

Besides the scanning radius pq, the formation length If is a second important parame- 
ter of CT physics. Changing energy u, one can vary If from If <^ Ra (quasi-instantaneous 
formation of the final-state hadron) to > Ra (the frozen-size limit) , and thus study the 
dynamical evolution of the small-sized, perturbative qq pair to the full-sized nonpertur- 
bative hadron. We predict very rich pattern of the nuclear shadowing and antishadowing 
phenomena, which changes with the scanning radius pq, i.e., with Q^, and with the evo- 
lution rate, i.e., with energy u. 

It is worth while to emphasize that the virtual photoproduction exemplifies in a partic- 
ularly transparent way the principle idea of CT tests of QCD: (1) A small transverse-size 
component of the interacting hadrons, p ^ Rh, is selected by the interaction dynamics. 
(2) The interaction cross section a{p) of the small-sized ejectile is measured by a strength 
of the final state interaction (FSl) in the target nucleus [4,5]. Notice a combination of 
the perturbative and nonperturbative aspects of QCD: the perturbative production of the 
small-sized ejectile is followed by probing its size in the nonperturbative, diffractive small- 
angle scattering in the nuclear matter. Nevetheless, QCD as a theory of strong interactions 
predicts that the strength of this diffractive scattering vanishes as p — >^ 0: cr (p) cx p^ [6,7]. 

The scanning mechanism. The quantum-mechanical description of the scanning 
goes as follows: At If, Ic > Ra the amplitude of the forward photoproduction on the 
free nucleon M^r and the nuclear transmission coefficient or the nuclear transparency 
Tva — daA/AdaN in the quasielastic photoproduction 7*^4 VA read [8,9]: 

Mr, ^ {V\a{p)\^*) (3) 

TrA^jJ d%m iyl^p^\,.^. (4) 

where T{h) = J dzuAip, z) is the optical thickness of a nucleus, the nuclear density n^(6, z) 
is normalized to the nuclear mass number A: / d^ruAif) — A. 
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The wave function I7*) of the qq fluctuations of the virtual photons was calculated in [3] 
in the mixed (p, z) representation, where 2; is a fraction of photon's light-cone momentum, 
carried by the quark, < 2; < 1. The most important feature of I7*) is an exponential 
decrease at large distances [3]: 

I7*) oc exp(-£p) (5) 

where 

s'' = ml + z{l - z)Q^ (6) 

Calculation of the matrix elements in (3,4) involves the (Ppdz integration. In the nonrel- 
ativistic quarkonium 2; pa 1/2, so that the relevant qq fluctuations have a size 

P ^ Po^ I = ~ , (7) 

What enters (3,4) is a product cr(p)|7*). Because of CT property, a{p) oc at small 
p [6,7], this product will be sharply peaked at p ~ p2 = 2pQ with the width Ap ~ 2pQ, 
which leads naturally to the idea of scanning: The transition matrix elements (3,4) probe 
the wave function of vector mesons at p ~ 2pQ, and varying pg by changing Q^, one 
can scan the wave function \V) from large to small distances. In Fig.l we demonstrate 
qualitatively how the scanning works. We also show the 2;-integrated wave functions of 
the ground state \V) and of the radial excitation \V'). 

The nuclear matrix element in (4) can be expanded in the moments (y\a{p)^\^*) which 
are the true QCD observables of the virtual photoproduction process. These moments 
probe the wave function \V) at different values of p ~ p2n- To the leading order in the 
final state interaction, 

TrA = l- ^vj I d%T{hf (8) 

where 

This expansion works well when 1 — Tta <^ 1, and is convenient to explain how the 
scanning proceeds. Notice, that a{py\^*) peaks at p ~ p4 = 4pQ. 



Scanning the vector mesons: the node effect. Start with the real photopro- 
duction (Q^ = 0) of the ground-state vector mesons (T, J/^', p°,...) and take the case 
of the charmonium. In this case pg = l/rric <^ Rj/^^ P2 is still smaller than i? ,//<!», but 
Pi ~ Rj/^i- For this numerical reason, one finds Sj/^ ^ atotiJ/"^ N), i.e., the predicted 
nuclear shadowing will be marginally similar [8,9] to the Vector Dominance Model (VDM) 
prediction [10]. This holds to much extent for the T and the light vector mesons as well. 
At larger Q^, one has p2, P4 <^ Ry and Ey ~ (t{pq) oc Pg with the calculable logarithmic 
corrections [3,11], so that the above marginal similarity with the VDM disappears. The 
nuclear transparency will tend to unity from below: 

l-Tr^oc-^pl (10) 
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The case of the radial excitation V is more interesting. Radial excitations have larger 
radius and larger free nucleon cross section: e.g., Utoti^'N) Ki 2.5atot{J/^ N) [8], which 
classically would suggest much stronger final state interaction for ^' than for J/\&. Pres- 
ence of the node in the V wave function leads to a rather complex pattern of the shadowing 
and antishadowing. In the photoproduction limit, because of p2 ~ -R-i-', there are rather 
strong cancellations between the contributions to the amplitude (3) from p below and 
above the node (the node effect, see Fig.l). For this reason, in photoproduction on the 
free nucleons, one predicts the ratio of the forward production differential cross sections 
r(g2 = 0) = da{-fN V'N)/da{jN VN)\t=o < 1- For the charmonium the ^'/(J/^) 
ratio r(0) = 0.17. In the T' production, the initial size pg scales as l/rUg, whereas the 
radius of the bb bound states decreases with niq less rapidly. For this reason for the T' 
the node effect is much weaker and we find the T'/T ratio r(0) = 0.84. The larger is Q^, 
the smaller size p2 is scanned, a contribution from the region above the node becomes 
negligible, and r{Q'^) will increase with up to r[Q'^ ^ my) ~ 1 (the exact limiting 
ratio depends on the wave function's at the origin). 

Since p^ is closer to the node position, the node effect is still stronger in the ma- 
trix element {V'\a{pY\-f). For the one finds (*'|cr(p)^l7) < 0' so that E*/ < and 
TrA{^') > 1 despite the larger frcc-nuclcon cross section [8,9,11]. In the T' production 
the scanning radius compared to the bottonium radius is relatively smaller than in the 
charmonium case, the node effect is weaker and FSI produces the shadowing of the T'. 
Still, in spite of atotC^'N) 3> atotC^N), shadowing of the T' is weaker than shadowing of 
theT. 

The (5^-dependent scanning changes the node effect significantly. The smaller the 
scanning radius pQ, the weaker are the cancellations. In the ^' case {V'\a{py\'j*) becomes 
positive defined, so that the antishadowing changes to the shadowing, which first rises 
with (5^, then saturates and is followed by the onset of asymptotic decrease (10). In the T' 
production the node effect is weaker, still it affects the Q^-dependence of the shadowing 
making it different from that of the T. 

Predictions for the heavy quarkonium photoproduction are shown in Figs. 2, 3. Because 
of the small size of heavy quarkonium, the results are numerically reliable, as they are 
dominated by the perturbative QCD domain. 

The two scencirios of scanning the light vector mesons. The case of the hght 
vector mesons is particularly interesting, as similar (anti) shadowing effects occur in the 
energy and range accessible at SLAC and at the planned European Electron Facility 
(EEF). Here at moderate "^m^ the scanning radius is large, pg ~ Ry, so that numeri- 
cally the predicted shadowing of the p° , Fig. 4, is not very accurate, although the accuracy 
increases gradually with as the scanning radius decreases into the perturbative domain 
Pq <^ Rv- Nevetheless, since the wave function of the p^ does not have a node, we beleive 
to describe correctly a smooth transition from the real to virtual photoproduction. 

In the photoproduction limit we find a strong node effect and strong suppression of 
the photoproduction of the radial excitation p' on nucleons, by more than one order of 
magnitude compared to the p° production, which broadly agrees with the scanty experi- 
mental data [10,12]. As here pq ~ Ry, we can not give a reliable numerical estimate of 
how small the p'/ p° cross section ratio is. Nevertheless, we can describe the two possible 
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scenarios of scanning the p' wave function, experimental tests of which can shed hght on 
the spectroscopy and identification of the radial excitations of the hght vector mesons: 

(1) The undercompensated free-nucleon amplitude: {p'\a{p)\'y) > 0. 

In this scenario the p' case will be similar to the case, apart from the possibility of 
anomalously strong nuclear enhancement Tta 3> 1, which might show strong atomic 
number dependence. Indeed, because of the larger relevant values of p the cross section 
(t(p) is larger, and the attenuation factor exp[— |(t(p)T(6)] in the nuclear matrix element, 
eq.(4), will suppress the large p region, effectively decreasing the scanning radius pq and 
diminishing the node effect. Detailed description of the atomic number dependence will 
be presented elsewhere. 

The (5^-scanning too will follow the ^'-scenario: change from the antishadowing to 
the shadowing with increasing Q^, followed by the saturation and then decrease of the 
shadowing according to eq.(lO). The range of at which the major effects should 
occur corresponds to a change of the scanning radius pq by the factor ~ 2, i.e., to 
Q"^ ~ Aruq ~ m?p. 

(2) The overcompensated free-nucleon amplitude: {p'\a{p)\'^) < 0. 

This scenario is preferred in the crude oscillator model used in [3] and fits the pattern of the 
node effect becoming stronger for the lighter flavours. In this case p2 > Ry and the higher 
moments too will be negative valued: (p'|o"(p)"|7) < 0, so that in the photoproduction 
limit one starts with the conventional shadowing: Tta < 1- 

In the (5^-scanning process the striking effect is bringing the free-nucleon amplitude 
{V'\(j{p)\Y) to the exact compensation at certain moderate Q^, when the decreasing 
P2 intercepts Ry (strictly speaking, because of the relativistic corrections and different 
quark helicity states, the compensation is unlikely to be exact). As a result, one finds a 
spike in TrA, Fig-3. With the further increase of one enters the undercompensation 
regime, and the futher pattern of the Q^-scanning will be essentially the same as in the 
undercompensation or the ^' scenario. 

The above described Q^-dependcnt scanning of the wave function of light vector 
mesons ofers a unique possibility of identifying the radial excitation of the light vector 
mesons (for the detailed discussion of the spectroscopy of light vector mesons see [12]). 
The corresponding experiments could easily be performed at SLAC and EEF. 



Quantum evolution, coherence and energy dependence of FSI. 

If the coherence length > Ra, then amplitudes of production on different nuclcons 
add up coherently. At moderate energy, Ic < Ra, the production rates on different 
nucleons at the same impact paramater b add up incoherently. In the opposite limit of 
Ic > Ra amplitudes of production on different nucleons add up coherently, and nuclear 
affects are generally weaker [8,9,11]. 

For the heavy quarks we have a strong inequality If ^ Ic [13]. The same inequality 
holds for the light vector mesons at ^ rriy. Consequently, at moderate energy, when 
Ic < Ra, one has still a broad energy range in which // > Ra and the transverse size of 
the qq pair is still frozen. In this regime the nuclear transparency is given by a simple 
formula 

TvA-jJd hdzuAih, z) lv\ain)\.^)^ ^^^^ 
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where t{b,z) = dz'nA{b, z'). Notice, that compared to the high-energy hmit (4), here 
the attenuation effect in the nuclear matrix element is weaker. The energy dependence of 
the nuclear transparency is given by an approximate interpolation formula 

Tr-Aii^) ~ Tr(/e « Ra) + F,^(/€)'[Tr^(/c > Ra) - Tr^(/c « Ra)] (12) 

where Fch{i<) is the charge form factor of the target nucleus and k = 1/lc = {Q"^ + niy) / 2v . 
(For the general idea of derivation of (14) and successful description of the NMC data 
[14] of the photoproduction of the J/^ see Ref.9). 

lilf < Ra, the spatial expansion of the qq pair becomes important. In can be described 
in either the quark basis used above, or in the hadronic basis, which is nice demonstration 
of the quark-hadron deuality. 

Consider the leading term of the final state interaction in eq.(lO) in the hadronic 
basis, i.e., in terms of Gribov's inelastic shadowing [15]. Inserting a complete set of the 
intermediate states, one can write down 

{v\a{pr\^*) = j:{y\^ipm{vMpw) (13) 

i 

In the hadronic basis the antishadowing of the ^' comes from the desctructive interference 
of the direct, VDM-hke rescattering 

Y ^^f' ^ ^' (14) 

and the off-diagonal rescattering 

7* ^ J/* ^ Psi' (15) 

(there is a small contribution from other intermediate states too). The reason for the 
strong cancellation is that the 7 — > transition is weak compared to the 7 — > J/^* tran- 
sition, whereas the J/^ — > ^' transition has an amplitude of opposite sign, and smaller, 
that the ^' — > ^' elastic scattering amphtude. To the contrary, in the J/^ photoproduc- 
tion, both amplitudes in the off-diagonal transitions like 7 — ^ \E'' ^ J/\E' arc small, and 
this explains why one finds a marginal similarity to the VDM in the photoproduction. 
At larger the node effect is no longer effective in suppressing the 7* — > transition, 
the off-diagonal amplitudes become significant for the J/ \E' photoproduction too and one 
finds strong departure from the VDM for the J/'^ too. 

By the nature of CT experiments, at finite energy (j(p)^ in the l.h.s. of eq.(13) is 
the nonlocal operator: the qq state produced on one nucleon is probed by a second 
nucleon a distance Az ~ Ra apart. As a result, the diagonal and off-diagonal rescattering 
amplitudes acquire the relative phase [15,16] = A.z{ml — m\)/2v ~ ^^/h- Upon 

the integration over only the intermediate states |i) for which A(/7ji ^ 1, or 

\ml - m^y\ < 2v / Ra (16) 

will contribute to the r.h.s. of eq.(13), so that the strength of the final state interaction 
will change rapidly from the near-threshold energy of <^ Ra to a higher energy of 
// > Ra- We describe the energy-dependence in this region using the quark-basis path- 
integral technique suggested in [3]. 
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At If « Ra the amplitudes of transitions (16) and (17) do not interfere. How- 
ever, since they are of comparable magnitude, for the ^' the VDM prediction for the 
nuclear shadowing breaks down even near the production threshold: Tr^ is larger than 
the Glauber model prediction. For the J/^ a similar incoherent contribution is small at 
small Q^, rises with as described above, and the near-threshold value of Tta rises too. 
In the ^' case the dominant effect of the Q^-dependent scanning is that the 7* — > ^' tran- 
sition amphtude increases with relative to the 7* — > J/^^ amphtude, which enhances 
the diagonal (shadowing) rescattering contribution compared to the off-diagonal (anti- 
shadowing) rescattering contribution. As a result of this competition the near-threshold 
value of TrA^^') first decreases significantly with Q^, followed by the J/^'-like behaviour 
at larger Q"^ (Fig. 2). For the both J/^ and ^' the near threshold behaviour of Tta is 
partially due to the kinematical rise of the threshold energy with Q"^. In the bottonium 
the intial size pq compared to the bottonium radius is relatively smaller, than in the 
charmonium, and Tr 4(T') exhibits a monotonous and v dependence (Fig. 2). 

With the rising energy the destructive interference of transitions (14) and (15) leads 
to a rapid rise of the nuclear transparency Tta and the onset of the antishadowing of the 
^' in the photoproduction limit. At the larger Q'^ the smaller p are scanned, and similar 
rise of Tta with energy ends up in the shadowing region. 

In the photoproduction of the light vector mesons If ~ If. and the more elaborated 
technique like the effective diffraction operator technique developed in [17] is called upon. 
This will be a subject of the further investigation. The effect on the Q^-dependence is not 
significant, though, and the cited results for the photoproduction of the p°-mesons, Fig.4, 
were obtained still assuming 1^ > If. At large Q^, when for all vector mesons If ^ Ic, the 
nuclear transparency for J/^, T exhibits a similar Q^- and z/-dependence. 

Coherence effects and scaling law for FSI. According to eq.lO, at large the 
FSI effect scales as A^/^/Q'^. However, this scahng law, suggested in [18], is valid only at 
the asymptotic energy [16]. 

Indeed, at moderate energy, the coherence constraint (16) limits the mimbcr of the 
interfering states N^ff. To a crude approximation, in this case the effective attenuation 
will be controlled by not a small size pg of the initial qq pair, but rather by the least 
possible size pmin of the wave packet, which can be constructed on a truncated basis of 
Ngff conspiring states [17].( Above we already have encounterd a strong effect of the 
coherence on the energy- dependence of Tta at fixed Q^.) Only in the very high energy 
limit, when N^ff is not bounded from above, Pmm ~ Pq- The constraint (16) is less strin- 
gent for the heavy quarkonia, as the level splitting Am <C 2mq, and is more noticeable 
for the light vector mesons, as here the m^-splitting between higher excitations, which 
enters the coherence condition (16), rises rapidly with the mass. In Fig. 5 we present our 
results for the dependence of 1 — TrA{p^) at fixed energy i/: it is much weaker than 
oc l/((5^ -l- rriy), which is an appropriate variable in blew of eq.(7). At asymptotic energy 
the scaling law (11) works for all vector mesons. Fig. 5. A somewhat late onsent of the 
scaling law (11) for the p° can be understood in terms of the overcompensation scenario 
in the p^, p' photoproduction, as one needs a relatively larger to make the node effect 
negligible. 
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Conclusions: We have shown that QCD observables of CT experiments correspond to 
scanning the non-perturbative hadronic wave functions with the Q^-dependent scanning 
radius pQ. The strength of the final state interaction in the exclusive virtual photopro- 
duction of vector mesons is shown to depend strongly on the nodal structure of wave 
functions. At large energy and we predict the universal pattern of shadowing for all 
the photoproduced vector mesons. The node effect in the virtual photoproduction can be 
used to identify the radial excitation states. The predicted Q"^ and energy dependence of 
virtual photoproduction on nuclei can easily be tested at SLAC and planned European 
Electron Facility. The dedicated experiments on virtual photoproduction of vector mesons 
desrve special attention, since theoretical predictions of CT signal are much more reliable, 
than in {e,e'p) or {p,p'p) reactions (for the recent review see [19]). 
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Figure captions: 

Fig.l. - The qualitative pattern of the the Q^-dependent scanning of the wave functions of 
the ground state V and the radial excitation V of the vector meson. The scanning 
distributions a(p)\l'^. (p) shown by the solid and dashed curve have the scanning 
radii pg differing by a factor 3. All wave functions are in arbitrary units. 

Fig. 2. - The predicted and i/-dependence of the nuclear transparency in the virtual 
phtotoproduction of the heavy quarkonia. The qualitative pattern is the same from 
the hght to heavy nuclei. 

Fig. 3. - The predicted Q^-dependence of the nuclear transparency in the virtual photopro- 
duction of the ^' (the undercompensation scenario) and p' (the overcompensation 
scenario) mesons. 

Fig.4. - The predicted Q'^ and i/— dependence of the nuclear transparency in the virtual 
photoproduction of the p°-mesons. 

Fig. 5. - Test of the scaling law (10) in the virtual photoproduction of the J/\]/ and p^ 
mesons: 

The right box - the p° production at fixed energy. The dashed straight hne corre- 
sponds to the 1/(Q^ + m^) dependence. 

The left box - at the asymptotic energy, the dotted straight line corresponds to the 
1/(Q^ + mj^^) dependence. 
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